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*""' Abstract 

Multiple Mellin-Barnes integrals are often used for perturbative calculations in particle 

r""! physics. In this context, the evaluation of such objects may be performed through residues 

'"^ calculations which lead to their expression as multiple series in powers and logarithms 

I of the parameters involved in the problem under consideration. However, in most of the 

i-G cases, several series representations exist for a given integral. They converge in different 

C^ regions of values of the parameters, and it is not obvious to obtain them. For twofold 

S integrals we present a method which allows to derive straightforwardly and systematically: 

(a) different sets of poles which correspond to different convergent double series repre- 
,_^ sentations of a given integral, 

J^ (b) the regions of convergence of all these series (without an a priori full knowledge of 

_- 1 their general term), 

C^ (c) the general term of each series (this may be performed, if necessary, once the relevant 

^— ^ domain of convergence has been found) . 

L.^^ This systematic procedure is illustrated with some integrals which appear, among 

,__! others, in the calculation of the two-loop hexagon Wilson loop in A/" = 4 SYM theory. 

T-H Mellin-Barnes integrals of higher dimension are also considered. 

> 
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1 Introduction 

The Mellin-Barnes (MB) representation is an important tool of asymptotic analysis and it 
is also very useful in high energy physics. In the context of the perturbative calculations of 
particle physics, for instance, the evaluation of a multiple MB integral, after resolving the e 
singularities, may often be obtained as a decomposition in terms of some MB integrals of lower 
dimension (typically onefold, twofold and threefold integrals) that have to be expressed as 
(multiple) series. Each of these MB integrals of lower dimension, which in general depend on a 
number of parameters equal to its dimension, do have several convergent series representations. 
These different series converge in different regions and are, as a rule, analytic continuations 
of each other. 

In this paper, we mainly focus on the case of twofold MB integrals (see however section El) 
and, after presenting a general method allowing to derive these different series representations 
from a given twofold integral without any use of analytic continuation properties, we show 
that for a large class of MB integrals met in the perturbative calculations of particle physics 
(if not all), one may derive the regions of convergence of such series without the full knowledge 
of their general terirQ Therefore, it enables one, if necessary, to choose the relevant region of 
convergence before a possibly extensive calculation. 

Although the method allows to get the complete series (and thus exact results) in the 
computation of a given quantity, it may also be interesting if one only needs the first few 
terms in the expansions of this quantity for different values of the parameters (which of 
course belong to different regions of convergence). Indeed, one does not need to perform the 
resummation of the series corresponding to one region, followed by an analytic continuation, 
in order to obtain the result in another region. 

We illustrate the presentation of our method with some examples of integrals which appear, 
among others, in the evaluation of the two-loop hexagon Wilson loop in AA = 4 SYM. 

2 Onefold Mellin-Barnes integrals 

In this introductory section we would like to briefly summarize the main different situations 
that may be encountered when one is interested in the series representations of onefold Mellin- 
Barnes integrals. To fix ideas, we begin with two simple illustrative examples. 

The first one is the toy integral Z{0) corresponding to the zero-dimensional version of 
the vacuum-to- vacuum generating functional of A (j)^ theory, used in [T] for the exposition of a 
method allowing to obtain non-perturbative asymptotic improvements directly from divergent 
perturbative expansions: 

Z(0) = -L / # e-i<^'-^*'. (1) 

It was shown in [T] that the following Mellin-Barnes representation for Z{0) may be obtained: 

c+ioo 
c—ioo 

where c = Res g]0, j[, and that by closing the contour of integration to the left, by Cauchy 



^We will see that this property is in fact not restricted to the case of double series. 



theorem, one obtains the asymptotic expansion 
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which, from d'Alembert's ratio test, diverges for any value of A. 

In jT], the quantity of interest (and starting point of the non-perturbative asymptotic 
analysis) was the divergent expansion ([3|, therefore it was not mentioned that one may also 
close the contour of pi) to the right to get the series 
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which, this time, converges for any value of A. 

Thus, ([3| and d4| are two very different series representations of Z{0). 

Let us come to our second example, the QED contribution to the anomalous magnetic 
moment of a lepton a^""^' ^° ' of the diagram of Figure hi 




Figure 1: 2- loop QED vacuum polarisation contribution to 5 — 2 of a lepton L. 



Following [2] or [3] it is easy to obtain the corresponding Mellin-Barnes representation 
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where c g]0, 1[. By closing the contour of integration to the right, one gets the serieaj 
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which, from d'Alembert's and Raabe-Duhamel's ratio tests, converges in the region —k- > 1. 
Alternatively, closing the contour to the left, one gets the series 
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which converges in the region —^ < 1. 

In fact, the two simple examples (pi) and ([5| reflect the particular situations that may be 
encountered in the computation of the following more general onefold Mellin-Barnes integral: 

c— joo 

where c is chosen such that the straight line of the integration contour does not intersect the 
poles of the integrand and Oj, bj, Ck and dk are real numbers. 
Let us define the quantities 

m n m n 

A = y^ Qj — y^ Ck and a = >^ \aj\ — N^ \ck\- (9) 

j=l k=l j=l k=l 

The region where the integral ([8| converges is | argx| < f a, see for instance |5 (one may 
sometimes extend the region of convergence on its boundary as for instance when a = 0). 
Three cases are in order |5], |6]: 

• If A > 0, closing the integration contour to the left leads to a series representation of 
the integral ([8| which converges for any value of x, but closing the contour to the righlrl gives 
a divergent asymptotic expansion. 

• If A < 0, closing the contour to the right leads to a series representation which converges 
for any value of x, but closing the contour to the left gives a divergent asymptotic expansion. 

• If A = 0, closing the contour to the left and to the right gives two convergent series (if 
the ttj are all of the same sign, one series is terminating), the first one converging in some 
disk |x| < R whereas the other converges outside this disk. Moreover, if a > 0, each of the 
two series is an analytic continuation of the other. 

It is easy to see that ([2| has A = — 1 and q = 3, and that ([5| may be reduced to the form 
^ with A = and a = 4. 

Therefore, the two expressions mj and ([7]), which may be expressed in terms of special 
functions as followa^ (with r = rnj/m^): 

^,.. r.. ^ ,-^V [_!_, + ililMl _ =^(i:M) + VArcCoth (^ In. - "^ 

^ r>i Vvr/ 4 Ar 4r 2 ^ -' 6 

+ -r ln(r) r ' ArcCoth (\/r) Inr — r Inll ] In r + r Li2 I - I 

(10) 

^In some cases there might be only a finite number of poles or no pole at all in this part of the complex 
s-plane, so that Cauchy theorem does not lead to a useful asymptotic information. This is the signal of an 
exponentially suppressed asymptotic expansion and other mathematical techniques have to be used. 

*Other equivalent representations may be found in the literature. 
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where $(z,s,a) denotes the Lerch's function ^{z,s,a) = Xln^o To+nF ^°^ l"^! ^ ■'■ ^^-^"^ '^ "^ 
0, —1, —2, ..., are analytic continuations of one another (see Figure^ where the point r = 1 
corresponds to the case of identical leptons I = L). 
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Figure 2: (^) Gj^"''^' ^° ' as a function of r: in blue, eq. (10) and in red, eq. Pll- 



In this paper, we are interested in the generalization of the results of this section to 
Mellin-Barnes integrals of higher dimension. We will be more particularly concerned with the 
extension of the A = case which seems to be the main situation met in the perturbative 
calculations of particle physics and we will see that for an integral which fulfills this property, 
several convergent multiple series representations may be obtained. As a rule, these series are 
analytic continuations of each other and we will show that it is possible to get the regions of 
convergence of these multiple series before computing their full expressions. 

3 Twofold Mellin-Barnes integrals 

In this section we consider twofold Mellin-Barnes integrals depending on two (complex) pa- 
rameters and we present a method which allows to derive straightforwardly and systematically: 



• different sets of poles which correspond to different convergent (double) series represen- 
tations of a given integral. 



• the regions of convergence of all these series (without an a priori complete knowledge 
of their general term) , 

• the general term of each series (this last step may therefore be performed, if necessary, 
once the domain of convergence of interest has been found) . 

This systematic procedure will be illustrated with some integrals which appear in the 
intermediate calculations of the two-loop hexagon Wilson loop in A/" = 4 SYM theory (see 
[8]), and on the example of the scalar box integral with one external mass (in the latter case 
we will also show how one can resolve the e singularities from a graphical point of view with 
this method). 

We saw, in the case of the two onefold integrals explicitly calculated in section [2| that 
the integration contours were closed to the left or to the right, the poles where one had to 
compute residues in the s-plane being either on the negative or on the positive Re s axis. 
We also saw that when A = 0, by closing the contour to the left or to the right, one obtains 
two convergent series. The generalization at the twofold level of this particular case does not 
lead to only two series but in fact to several. These convergent (double) series are, as a rule, 
analytic continuation of each other and the poles where one has to take the residues which 
correspond to each series belong to some specific sectors in the (Re zi, Re Z2) plane, where z\ 
and Z2 are the two integration variables [9] . Of course in general these sectors (the so-called 
cones, in the following) are not accessible by naively closing the two contours of the twofold 
integral to the left or to the right, but we will show that it is anyhow easy to find the cones 
by simple intuitive rules or with the help of a geometrical procedure. 

To expose the method in all details, we will first consider the simple twofold Mellin-Barnes 
integral i?_i(ni, n2, Us) (first line of eq. (3.6) in [8]). We will show that for this integral five 
different cones may be obtained and, after recalling the way that residues which correspond 
to a given cone are calculated, that each of the five associated double series converges in a 
particular region of the (|mi|, |n2|) plana^that we will give explicitely. Then, we will detail 
how to get two of these five series from their respective cone. 

We insist to say that one of the important points of this paper is that the region of 
convergence of a given series may be obtained without the full knowledge of the general term 
of this series. Therefore, in the cases where one knows the numerical value of the expansion 
parameters (like it is, for instance, in the QED calculations done in [3], where the parameters 
are lepton mass ratios) this allows to select the relevant cone right from the beginning without 
having to compute everything. However, we will also see in the example of R-\{ui,U2) that 
there exists a (semi) infinite region in the (ui, U2) plane that is not accessible to any of the five 
series representations that our method allows to find. Similar "out of reach" regions will also 
be met in the case of the two other integrals under study in this section. The problem of the 
direct obtaining of a series representation which converges in these (at first sight) inaccessible 
regions (without using the traditional approach, i.e without performing analytic continuation 
of the resummation of the known result corresponding to one of the accessible regions) is 
presently under investigation with a complementary original method but will not be treated 
in this paper. 

After the calculation of i?_i(ui,U2) we will consider in this section the slightly more 
complicated twofold Mellin-Barnes integral i?_i_2i('Ui,U2)'U3) (second line of eq. (3.6) in |8]) 



^Indeed, in can be seen in [^ that 7?_i(ui, 112, us) does not depend on M3, therefore in the following we do 
not keep the notation of these authors and, instead, we will write i?„i(iti, 112). 



and the scalar box with one external mass. This will allow us to explain, among other things, 
a technical point of the method which may be important to underline for the reader and 
which concerns "spurious" sets of poles which, if misinterpreted, could be thought as being 
part of the set of relevant poles inside a given cone. These undesired singularities should 
however not be taken into consideration in the calculations and we will see that it is easy to 
distinguish them. Moreover we should add that if one were anyway trying to compute them 
by mistake, one would fall on undecidable ill-defined results which would therefore signal 
their presence. Apart from this fact, the i?_i_2j(tii,U2,«3) and scalar box integrals are not 
harder to compute than i?_i(ui, U2) since it is worth recalling that our method is a systematic 
one. In particular, we will not need any introduction of Euler integral representation or other 
integral manipulations as was performed in |l8J: keeping right from the beginning the original 
twofold Mellin-Barnes representations and applying our method directly on them allows to 
obtain the different convergent double series straightforwardly. 

3.1 The R_i(ui,U2) integral 

In the following we consider ui and U2 as complex numbers. Let us start by giving the 
definition of i?_i(ni,U2) (first line of eq. (3.6) in ^): 

R-l{ui,U2) 

c+ioo d+ioo 

= - / ^ / ^ <^ ^r r' (-^1) r (^1) r^ (-^2) r (z^) r (zi + ^2 + 1) , (12) 

c— zoo d—ico 

where c = — o and d = —^. Note that there is a symmetry of this integral by exchanging ui 



and U2 since we could take c = d = —h orc = d=— 4 without changing ( 12 ). 



It will of course allow us to obtain some of the series representations directly from others. 



The singular structure of the integrand of (12), in the (Re zi,Ke Z2) plane, is represented 



in Figure pi It is easy to see that, with n G M, it corresponds to a drawing of the (singular) 
lines —zi = —n (the vertical blue lines on figure ^ and —Z2 = —n (horizontal blue lines) 
which both appear twice, of the lines zi = —n (vertical red lines), Z2 = —n (horizontal red 
lines) and of the lines zi + Z2 + 1 = —n (oblique green lines). 

At the intersections of these singular lines are the singular points where one will have to 
compute residues of the integrand, the resummation of these in the different coneaj leading 



to various convergent (double) series representations of (12). 

It is clear from Figure p] that there is only a finite number of intersection types (and a 
finite number of domains in the (Re 2;i,Re Z2) plane where one will always find the same 
intersection type), which makes the resummations possible. Moeover, apart from a few points 
as for instance the point (0, 0), which is by itself its own type, to a given type of intersection 
correspond one series of residues. For instance, as we will see, in the highest right-corner of 
Figure p] for (zi, Z2) = (1 -|- m, 1 -|- n) where ?Ti > and n > are integers, the intersection 
type consists in a crossing of two horizontal and two vertical (blue) lines and corresponds to 
a double series. Still in this corner, for {zi,Z2) = (0, 1 -|- m) (resp. (^1,^2) = (1 + m,0)) 
where m > is an integer, the intersection type consists in a crossing of two horizontal 
blue lines (resp. two vertical blue lines) and three vertical lines, one red and two blue (resp. 
three horizontal lines, same colors). The latter correspond obviously to single series, since 



As we said above, there are cases where some of the intersections are irrelevant, but this does not occur 
in the computation of R-i. 
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Figure 3: Singular structure of the integrand of R^i{ui,U2)- Coinciding singular lines have 
been slightly separated for convenience. The red dot is the point (c, d). 

the parametrization of these singular points do only depend on the m summation indice. For 
(zi,Z2) = (0,0), the intersection type consists in a crossing of three horizontal and three 
vertical lines (one red and two blue in both cases). In fact, if one takes the residues at the 
intersections which belong to this complete quarter of plane, one matches exactly equation 



(C.2) of [8j which has been obtained by closing both contours of (12) to the right. 



We will show now that one can obtain this cone as a particular result of the application 
of a simple rule which allows more generally to exhibit various other sets of singular points 
from a given twofold Mellin-Barnes integral {i.e. other compatible cones), corresponding to 
other convergent series representations of the considered integral. 

3.1.1 Determination of the cones 



First, we have to find the fundamental polygone of (12), i.e the region of convergence of the 



integral in the (Re zi,R.e Z2) plane where a continuous change of the point {c,d) from its 
initial value (—3, —4) may be performed without modifying the value of the integral. It is 
dictated by the simultaneous constraints Re zi < 0, Re Z2 < and Re(l + zi + Z2) > 
(leading to the red triangle in Figure H|. 

Then, by a simple inspection of the constraints implied by the different coupleqjof gamma 
functions of the integrand of (12), whose singular lines intersect each other in some regions 



of the (Re zi,Re Z2) plane, it is straightforward to get the different cones. Indeed, from the 
eight different couples of gamma functions which may be obtained, only five lead to non-trivial 
regions which are compatible with the fundamental triangle, in the following way: 

• (r (— ^i) , r {—Z2)) leads to Cone 1: Re zi > 0, Re Z2 > 0. 

• (r (zi) ,r (— Z2)) leads to Cone 2: Re zi < —1, Re Z2 > 0. 



^A given couple of gamma functions should of course not depend on only one of the variables zi or Z2. 



• (r (2:2) ,T{zi + Z2 + 1)) leads to Cone 3: Re Z2 < -1, Re(l + zi + Z2) < 0. 

• (r {—zi) , r {Z2)) leads to Cone 4: Re zi > 0, Re Z2 < — 1- 

• (r {zi) ,r{zi + Z2 + 1)) leads to Cone 5: Re zi < -1, Re(l + zi + Z2) < 0. 

The first three cones are depicted on the first picture of Figure [4j the two others on the 
second picture. 
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Figure 4: In grey, the five different cones of i?_i(ni, U2)- 

However, it may be more convenient for the reader to know that the cones may also be 
obtained following a simple general geometrical procedure |9|. As described in this reference, 
one has to consider the value of the vector 






(13) 



which generalizes the quantity A of section [2] and where the real coefficients Oj , bj , c^ and 
dk are defined in the general twofold Mellin-Barnes integral as follows 
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c+ioo d+ioo 
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where c and d are chosen such that the integration contours do not intersect the poles of the 
integrand. 

Now the statements at the end of section p] may be generalized as follows |9] : 

• If A = (0, 0), we are in a so-called degenerate case where several relevant cones coexist, 
corresponding to different convergence domains. Therefore, several convergent double series 
expansions may be obtained which, as a rule, are analytic continuations of each other. 



• If A 7^ (0,0) and none of the vectors {aj,bj) is proportionarl to A, then we are in a 
simpler non-degenerate situation where there is only one double series representation which 



*The case where A 7^ (0,0) and at least one of the {aj,bj) is proportional to A is also a degenerate case 
(in [To] a quantum field theory calculation with this particular type of degeneracy has been considered). 



converges absolutely for all non-zero complex values of x and y. 



Looking at (12), one concludes that i?_i(ni,n2) has A = (0,0). In fact, it seems that 
in the high energy physics perturbative computations, one often faces twofold Mellin-Barnes 
integrals with A = (0, 0) (in this paper we only consider this case). 

Let us now come back to the geometrical procedure that we mentioned above for the 
determination of the different cones, in order to describe it in detail. For this one has to draw 
on Figure Is] an arbitrary straight line I that goes through the point (c, d), in our case the point 
(— Tj, — j). Obviously there is an infinity of such lines, but only a finite number of compatible 
cones can be deduced from this infinite set. 

Once the straight line / is drawn, the extraction of the corresponding compatible cone(s)r| 
follows from the following rule: the singular lines that cross I on one of its "sides" (the point 
{c,d) defines the separation point between each of the two "sides" of /) intersect the singular 
lines that cross I on its other side at some points which are inside the compatible cone [9]. 

In Figure [5] we show one possible straight line from which two compatible cones emerge 
(the colored regions), the right 90 degree cone corresponding, as we already said, to the results 
given in eqs (C.2), (C.3), (C.4) and (C.5) in Appendix C of [8]. 
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Figure 5: Cone 1 and Cone 5. 

Straightforwardly, the three other cones of Figure |4] may be obtained by a rotation of the 
^-line around (c, d). 

All cones being found, the integral i?_i(ni, U2) may then be computed as follows (we use 
some of the notations of |9] for later purpose) : 
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^In some cases one may find two such cones for a given I: one in each of the two half-planes that are 
separated by I. 



where w is the meromorphic 2-forni 

w = u{^ ui^ r^ {-zi) r [zi) r^ (-^2) r (^2) r (zi + Z2 + 1) dzi a dz2 , (i7) 

7 = (c, d), and where Res [.] is meant for residue. 



The sum in ( 16 ) has to be understood as a sum over the residues taken at the intersections 
inside the i^^ cone and TZi is the region of convergence of such a sum. 

Let us now explain how to get the different regions of convergence TZi associated to each 
of the five cones. 

3.1.2 Regions of convergence 



The different cones corresponding to ( 12 ) have been found and an obvious need now is to know 
in which region of the (|ni|, |m2|) plane the corresponding double series converge. Although 
it is possible to get the region of convergence of a given series once one has its general term, 
this order of doing things is not always interesting since in many instances one knows right 
from the beginning the numerical values of the expansion parameters, and therefore one would 
prefer to know to which region of convergence these values belong, in order not to compute all 
the cones one by one but only the "right-one" directly. We will show that it is possible to find 
the different regions of convergence without a full computation of the general terms of the 
corresponding double series, thanks to a nice property of the integrand of the Mellin-Barnes 
integrals met in high energy physics perturbative calculations. To explain this in full details, 
it is now necessary to make a digression in order to define several quantities that will be 
needed. 

Residues computational method Inside a given cone, the first important thing to do is 
to consider exclusively the intersections where at least one singular line crosses one side of 
the line / and one singular line crosses the other side of I (the sides of / have been defined 
in section 3.1.1 ). Indeed, it may happen that, for a given cone, part of the intersections of a 
singular structure are formed only by singular lines which cross the same side of I. These are 
the spurious singularities that we talked about in the introduction of section [3] The reason 
why these intersections should not be considered is simply that they do not comply with the 



definition that we gave in section 3.1.1 for a singularity to belong to a cone, and we prefer to 
insist on this point in order to prevent the reader from an incorrect use of the method. In the 
calculation of R^i(ui,U2) no such spurious singularities will be met, but as we already said 
this will happen for some of the cones of R-i-z^ (ui, U2, its) (the small black circles in Figure 



10 distinguish these spurious singularities) and in many other instances as for example the 



scalar box integral discussed later in this section (see the small black filled circles in Figure 



13) or the QED calculations that have been performed in fj]. 

Once the possible spurious singularities have been discarded from the calculation, one has 
to classify the different subsets of relevant singularities (indeed, we have seen that inside a 
given cone, there exists in general several types of intersections). This way, one may deal 
with infinite subsets of residues at the same time, which may be convenient if the aim is to 
get an exact result. If only the first few terms of the expansion are needed, they come from 
the contributions of the intersections which are the closest to the fundamental polygon and, 
of course, they may be computed one by one (this was the strategy adopted in [3]). 

The coordinates of a particular subset are parametrized by affine functions / and g as 

{zi, Z2) = {f{m, n),g{m, n)), (18) 

10 



where m and n are integers. / and g are trivially obtained either by looking at the picture 
representing the singular structure of the integrand (Figure^ in our present case) or by solving 
the system of the singular lines equations coming from two gamma functions chosen between 
all the gamma functions which are involved in the considered type of intersections. 

Then, for the computation of the corresponding residues, the following four steps have to 
be performed [3l [11] : 



One first performs on (17) the change of variables zi i— ?■ zi+f{m,n) and Z2 ^ Z2+g{m,n) 



so as to bring the singularity to the origin (0, 0). 

• One then uses the generalized reflexion formul"'^'' 



n-k + z) = (-1)^ -^ ^ = (-1)^ r(i + z)r(i-z) 

(where /c > is an integer) in order to explicitly extract, in the denominator, the singular 
structure of the singular gamma functions. 

• Two different cases may then appear, depending on the type of intersection considered 
(or, in other words, the type of arguments of the singular gamma functions which contribute 
to this kind of singularity) : 

First case: if the singular lines crossing each other at the intersections are only horizontal 
and vertical lines (i.e if the arguments of all the singular gamma functions do not depend on 
zi and Z2 simultaneously) , one finds for the 2-form an expression of the type 

^='^^dz,Adz2, (20) 

where n and m are positive integers and the function h is analytic at the origin. 

Second case: if at least one of the singular lines is neither vertical nor horizontal or, 
equivalently, if there are some oblique lines at the intersection points (z.e if the argument of 
at least one of the singular gamma functions depends on zi and Z2 simultaneously), one finds 
the more general expression 

W = ^— 5/''^'J ■,, r s dzi A dZ2 , (21) 

[zi + a Z2)P[zi - b Z2P z\ zl 



where a, h are positive numbera^ and p, g, r and s positive integers (possibly null for q if 
p 7^ or conversely, idem for r and s). 



To find the expression ot the residue of the form (20) one may directly apply the Cauchy 
formula 



D . 1 1 a"+'"-2/l(zi,Z2) 

Res [wj — 



(n-l)!(m-l)! dz^-^dz 



ni~\ 
2 



n — 1 ! 771 — 1 ! 

(0,0) ^ ' ^ ' 



'^''One may also use the identity r(— ti -^ z) — ^j-j„' t^L^-. but it is often more convenient to use (|l9|. 

^^ Notice that in the case of the MeUin-Barnes integrals met in perturbative quantum field theory these 
positive numbers are in fact integers. It is important since this is a condition required by Horn's theorem to 
find the regions of convergence of the series (see Appendix p| . 
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but one has to care about the fact that an overall sign may contribute, related to the (clockwise 
or anti-clockwise) integration contours. As we will see, to find this sign it is more convenient 



to consider (20) as a trivial case of the more general situation (21 ) witli p = q = 0. 



On (21) one cannot trivially apply the Cauchy formula, one first has to use the so-called 



transformation law |121 l3l ITT], The transformation law applied to (21) reads 

h{zi, Z2) detA 



Res 



h{zi,Z2) 

(zi + az2)P{zi - bz2y z\ z. 



dz\ A dz2 



Res 



dz\ A dz2 



(23) 



if one can find a 2 x 2 matrix A whose elements are analytic and which satisfies Aj = g 



where g 



{Zi,Z2J 



and j are positive chosen integers) and / = (/i,/2) , /i being the 



product of a speciyic subset ^^ of the following terms: {zi-'raz2)'^, {zi—bz2y, 2;[ and z|, while 
/2 is composed of the product of the terms which do not belong to /i (in some cases one of 
these specific subsets may be reduced to only one term, then there is of course no need for a 
product). Moreover one has to check that /~^(0,0) = (0,0) and g~^(0,0) = (0,0). Once this 



has been done, (22) may be applied on the right hand side of (23). 



In fact the expression of /i and /2 are precisely defined in Figure |6] which summarizes the 
two different possibilities. One can see that the red direction vector defines the orientation of 



Re Z2 



Re 22 



Re zi 





Figure 6: Definition of /. 

the Hine which allows to consider two half-planes: the "West" and the "East" half-planes (the 
direction vector gives the "North"). If the cone under study belongs to the East half-plane, 
all singular lines which cross the /-line on the side where the arrow of the direction vector lies 
contribute to f2- All others singular lines contribute to /i. If the cone considered belongs to 
the West half-plane, this is the converse situation. 

The correct sign of the residue, which is automatically given with this approach in the case 



of (21 ), is then also easily obtained for Cauchy type residues (20 ) by the value of det A = ±1 



in the trivial application of the transformation law 

1 0^ 
.0 1, 



(24) 



We stress that this sign is undecidable for spurious singularities of the type ( 20 1 when considered as a 



trivial case of (21 1. This reflects their spuriousity. 
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This point is very important since an arbitrary subset gives wrong results. 
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or 



1 

1 



(25) 



depending on the vector / that has been found. 

This definition of the transformation law explains why one cannot decide for the value 
of the residues in the case of the spurious singularities that we mentioned above: in this 
particular case the / vector cannot be obtained since the singular lines cross the Hine only 
on one of its side. Therefore the transformation law cannot be applied, even in its trivial 
pl) or (pSl). 



versions 



In our present case of study, it is clear that for all residues associated to the upper cone of 
Figure p] (Cone 1), a direct application of Cauchy formula is possible, the sign of the residues 
being fixed by eq. (24). On the contrary, in the lower cone (Cone 5), one will have to apply 
a "non-trivial" transformation law as an intermediate step in the calculation of each of the 
different contributions. 



• At the end, the partial derivatives of (22| or (23 1 may be easily obtained by hand (see 
Appendix A for a formulary and some notations used in the calculations) in terms, among 
others, of polygamma functions tp^^' of different orders which may be translated as harmonic 
numbers H]\j = X^j^^ ^ and their generalizations H^^ = X]i=i ^ |13| if desired, using the 
relations 

i;{N + 1) = H{N) - 7E (26) 



and 



^ W (AT +!) = (_ l)'=+iA;! 



CfcH 



H 



(fe+i) 

N 



(27) 



where A^ and k are integers. 

The infinite double sums of the final results may in turn be expressed in terms of nested 
harmonic sums if one wants to rewrite them as a combination of multiple polylogarithms with 
the help of the different algorithms developed in [E 



Back to the regions of convergence Remember that our aim is to find the region of 
convergence associated to a given cone without a full calculation of the corresponding double 
series. For this, we need some of the tools described in the last paragraph. 

It is easy to see that from the use of (19), after the change of variables performed on w, 
the h function in (20) or (21) has the form (130) with x=—,y=—,s = zi,t = Z2 and 



where ej and fk arehnear combinations of m and/or n. Moreover, from (22) as well as eqs 



(132), (137), (138), (139), (140), ( |141 ) and their higher-order derivatives, one can infer that 
any residue may be written as /i(0, 0) multiplied by a combination of polygamma functions 
of different orders. The latter depend also on the summation indices m and/or n. The case 
where both m and n are involved in the argument of the gamma functions of /i(0, 0) is in 
general the most restrictive concerning the regions of convergence and corresponds of course 
to the occurence of a double series. We proved in Appendix [B] section B.4| that the polygamma 
functions do not affect the region of convergence that would be obtained by simply considering 
that the general term of the double series is /i(0, 0). 

Therefore, since h{0, 0) is a ratio of gamma functions, and since in the quantum field 
theory perturbative calculations the coefficients of the integration variables in the argument 
of these gamma functions are integers, one may use Horn's theorem to find the regions of 



convergence (see Appendix B, section B.2.1 ). We show how to proceed in the following. 



13 



Cone 1 To Cone 1, as we said above, corresponds only one double series, coming from 
the residues of the poles at (1 + m, 1 + n). We then have 



/i(0,0)| 



(l+m,l+n) 



u'^+\^+^ r(3 + m + n) 
(1 + m)(l + n)r(m + 2)r(n + 2) 



By eqs. (142) and (143) we find 



and 



F(m, n) 



G(m, n) 



m 



n 



m 

m + n 



n 



(28) 



(29) 



(30) 



Since we are in the p = r and q = s case (in the notation of Appendix p] section B.2.1), the 
region of convergence of the double series is then given by TZi = C (I'D where the domains C 
and V are defined as (see Appendix IbI) 



C= \ (|'Ui|,|u2|) I < |ni| < 



1 



|F(1, 



and < 1^2! < 



1 



|G(0,1)| 



and 



V= \{\ui\,\u2\) I V(m,n) G R^ : < |ni| < — ^ or 0< jusl < ^l 

[ ' |F(m,n)| |Cj(m,n)|J 



Therefore, one has 



and 



C = {0 < \ui\ < 1 and < \u2\ < 1} 



(31) 



(32) 



(33) 



V = {\ui\ + \U2\ < 1}. 



Indeed, (34) fulfills the disjunction constraints of (|32|): if < \ui\ < yp? m 



(34) 

and < \u2\ < 
.„, — yj it is straightforward to prove that (29) and (|30) lead to |ni| + \u2\ < 1- Now suppose 
;hat < luil < t^, ^t is false, i.e that lui > j^, y,. Then, inserting this in the relation 

I -^1 |F(m,,?i)| ' I ^ |F(m,n)| ' <= 

^i| + 1^2! < 1 one obtains, still by (29) and (^30), that < \u2\ < ^ 



|G(m,n)|- 



The single series corresponding to the poles at (1 + m, 0) and (0, 1 + m) do not lead to 
more restrictive constraints since we have 



/i(0,0)|(i+„o) 



u 



m+l 
1 

1 + 771 



and 



MO, 0)1 



rn+l 



(0,1+m) 



(35) 



(36) 



1 + m 

from which, by d'Alembert ratio test, one obtains respectively {\ui\ < 1} and {\u2\ < 1}- 

The region of convergence of Cone 1 is then given by TZi = {\ui\ + \u2\ < 1}, and it is 
represented in yellow in Figure [7j Let us remind that we did not need the full computation of 
the general term of the double series, which is given in (54 ), to get its region of convergence. 

From the symmetry of the integral i?_ 1(1*1,112) that we mentioned above and which is 
also visible in Figure |4j it is obvious to deduce that the region of convergence of Cone 2 may 
be obtained from the one of Cone 4, and that the region of convergence of Cone 3 may be 
obtained from the one of Cone 5 (or conversely). 



14 



Ilial 1-5 




7^, 



0.0 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 7: Regions of convergence of the series representations of R-i{ui, ^2). 

Cone 4 and Cone 2 In the case of Cone 4 we have three types of intersections lying 
at (1 + m + n, —1 — m), (1 + m, —2 — m — n) and (0, —1 — m). Therefore we have two double 
series and one single series. 

The single series gives the constraint {\u2\ > 1}. The series associated to the poles at 
(1 + 771 + n, — 1 — m) has 



^(0)0)l(l+m+n-l-m) 

Therefore 

and 

and we find 
and 



i-iy 



m+l 



r(l + m + n)r2(l + m)r(l + n) 
r2(2 + m + n)r(2 + m) 



m 



m + n 
n 



G(m,n) 

m + n 

Ci = {0 < |ni| < 1^2! and < \ui\ < 1} 



(37) 

(38) 

(39) 

(40) 
(41) 



Vi = {\ui\-1< \U2\}. 

Using the same procedure, one finds for the last series (poles at (1 + m, —2 — m — n)) 

C2 = {0 < |ni| < \u2\ and 1 < |n2|} (42) 

and 

V2 = {\UI\+1<\U2\]. (43) 

The region of convergence of Cone 4 is then 7^4 = Ci Pi P2 (the orange region in Figure u\ . 
This immediately gives the region of convergence of Cone 2 (the green region in Figure u\ 



as 



7^2 = {0 < \u2\ < \ui\ and < \u2\ < 1} n {\u2\ + 1 < \ui\} . 



(44) 
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Cone 3 and Cone 5 The same procedure applied to Cone 3 gives 

7^3 = {1 < \ui\ and 1 < |m2|} n {|mi| + 1 < |n2|} (45) 

(the cyan region in Figure W\ from what we deduce that 

7^5 = {1 < \ui\ and 1 < |n2|} n {|n2| + 1 < |ni|} (46) 

(the magenta region in Figure Q. 



In conclusion we see that the symmetry of the integral ( 12 ) is of course also reflected in the 
regions of convergence of Figure [7] However one may notice that a semi-infinite white band 
is not reachable by any of the series representations that our method allows to obtain. In fact 



we believe that it is possible to extract, from the Mellin-Barnes representation (12), other 
series that would probably fill this gap, from a complementary method which is presently 
under investigation. 

We checked that the regions of convergence TZi, ..., TZ^ may be obtained from the results 



of section B.5 since the corresponding series (with general term /i(0, 0)) are in fact Kampe 
de Feriet double series. 

Let us now apply the computational strategy described above to obtain the different series 
representations which converge in each of the regions of Figure [7] 

3.1.3 Cone 1: series representation valid in the region TZi = {\ui\ + \u2\ < 1} 



We begin with the simplest case, namely Cone 1. As we said in section [37L| there are four 
types of intersections that have to be considered separately and whose contributions will be 
subsequently added. From the discussion on the transformation law, we also concluded that 
there is no need for "non-trivial" applications of the transformation law in this cone and, 
moreover, that no additional overall sign appears for every residues associated to this cone 
since Figure |6] implies that for Cone 1 the sign is fixed by the determinant of the matrix in 



(24). 



Type 1: the singular point (^1,^2) = (0,0). 

In this case obviously no change of variable is necessary and one therefore has 



^ I Type 1 = ^3^3 ^^1 ^ dZ2 (47) 



I Cone 1 hi{zi,Z2) 

where 

hiizi, Z2) = -ul^ul^ r^ (1 - zi) r (1 + zi) r^ (i - ^2) r (1 + Z2) r (zi + ^2 + 1) • (48) 

Using the Cauchy Formula ( |22[ ) one obtains 

„ / N I Cone 1 

= ^/iP)(0,0) 

= - f -In^ui +C2 j (In2u2 + 4C2) - C2lnnilnn2 -C3(lnni -Mnu.2) - -C| - -C 
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4 • 



(49) 



Type 2: the singular points (zi, Z2) = {1 + m,l + n) where m,n > are integers. 
Performing the change of variables zi h^ zi + 1 + m, Z2 >-^ Z2 + 1 + n one has 



CO 



I Cone 1 
Type 2 



-U 



zi+l+m 22+l+n ■p2 



r^ (-Z1 - 1 - m) r (zi + 1 + m) 
X r^ {-Z2 - I - n)r {z2 + 1 + n)r {zi + Z2 + ^ + m + n) dzi A dz2. (50) 



After applying (19), one obtains 



UJ\ 



I Cone 1 h2{zi,Z2) 



Type 2 



2^2 



ztz, 



dzi A dz2 



(51) 



1^2 



where 



h2{zi,Z2) = -ul-^^+^ ^.^^+1+" ^' ^^ + ^^^ ^' ^^ " ^^^ 



Uo 



2 T^{zi+m + 2 

T\l + Z2)T\l-Z2) 



r (zi + 1 + m) 



r2 {z2 + n + 2) 
Using now the Cauchy Formula one has 



T{z2 + l + n)T{zi + Z2 + 2, + m + n) 



(52) 



00 00 

i?_i(ni,n2)i^:::j = ^j:4''')(o,o) 



I Type 2 



(53) 



m=0 n=0 



and from equations (137), (133), (134) and (135) or (136) a straightforward (mental) calcu- 
lation gives 



R^l{ui,U2) 



Cone 1 



I Type 2 
°o 00 ,m+l„,n+l 



^ (l + m)(l 



'u™+^M^+^ r(3 + m + n) 



m=0 n=0 



)(l + n)r(m + 2)r(n + 2) 



In ui + ^(1 + m) + "0(3 + m + n) — 2-0(2 + m) 



In ^2 + V'(l + "-) + "^(3 + m + n) - 2-0(2 + n) 
+ V'^^^(3 + m + n)} . (54) 

It may be useful to express the polygamma functions in terms of harmonic numbers and their 



generalizations using the relations (26) and (27). 
One then finds 



Cone 1 
Type 2 
°° 00 m+l,,n+l 



R-l{ui,U2 

00 c 



m=0 n=0 



v^+'u^-^' r(3 + m + w) 
(1 + m)(l + n)r(m + 2)r(n + 2) 



Hm + H2+m+n 



n + 1 



m + 1 

- -f^n + H2+m+n + In ^2 



+ Inui 



-I- 1,2 -n2+m+„ 



(55) 
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Type 3: the singular points (zi, Z2) = {1 + m, 0) where ?ti > is an integer. 



Performing the same change of variable as in Type 2 for zi and no change of variable for 
Z2 one has 

iConel _ h3{zi,Z2) 



UJ 



Type 3 



2^3 



ZtZ, 



dzi A dz2 



(56) 



1^2 



where 



^3(^1,2:2) = -ul 



21- 



+1+™^,, r2(i-zi)r2(i + zi^ 



u 



r(zi + l + m) 



r2(2 + m + zi) 

X r^ (1 - Z2) r (1 + Z2) T{zi + Z2 + 2 + m) , (57) 



and, by the Cauchy Formula (22) and equation (141), one obtains 



n / \\ Cone 1 

i?_i(ui,n2)|^yp^3 



^ 00 
m=0 

^ 00 ?Tl+l 

y^ ^ /[Inni + V(l + m) - ^(2 + m)] ([Inns - V'(l) + ^^(2 + m)f + 3V'^^^(1) 

2 ■'^— ^ 1 + 771 •- V 

+^^^^^(2 + m)) + 2 [lnn2 - V(l) + "^(2 + m)] ij^^\2 + m) + ^^^^(2 + m)\ 



m=i) 



m=0 V L \ / 

+ (inui - ^) (4C2 - H^X) - 2 (C3 - ^Si) } . 



[ In U2 + -f^m+l] 



(58) 



Type 4: the singular points (21,^2) = (0, 1 + rn) where rn, > an integer. 

Because of the symmetry of the integral R-\ this case may be obtained from Type 3 by 



exchanging u\ and U2 in (58). 



As a numerical check of our results, we directly compute the integral (12) for u\ = 0.1 
and U2 = 0.2 and find 

i2_i(0.1, 0.2) ~ -37.976338. (59) 

Now, a truncation of the series expansion of Type 2, 3, 4 at M = A^ = 10 added to the 
contribution of Type 1 for the same values of the arguments gives the following result 



I Cone 1 



I Cone 1 



^-l(0-l'0-2)|T;;:2,M=;V=10 + «-l(0-l'0-2),T,pe3,M=10 

+^-i(o-1'0-2)i?:::;m=io+«-i(o-1'0-2)It:::i 

~ -0.3035221 + 1.0827389 + 2.2579423 - 41.0134972 
~ -37.976338. 



(60) 
(61) 



We also checked that eqs. (C.3), (C.4) and (C.5) in Appendix C of [8j may be obtained from 



eqs. (49), (55) and (58) and its symmetrical expression. 
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3.1.4 Cone 3: series representation valid in the region 7^3 

{\ui\ + 1 < 1^21} 



{1 < \ui\ and 1 < |n2|}n 



We now turn to a less trivial cone, in the sense that the intersection types all involve oblique 
lines and therefore imply a use of the transformation law. This adds to the computation an 
intermediate step, but apart from this the calculation of Cone 3 follows exactly the same 
reasoning as for Cone 1 and the complexity is not increased. 

Let us begin by looking for the different types of intersections in Cone 3. It is clear from 
Figure [8] that there are 3 different types. 

Rez2 




Rezi 



Figure 8: Cone 3. 



Type 1: the singular points (21,22) = ( — 1 



m, 



-1 — n) where m,n >0 are integers. 



has 



Performing the change of variables 21 ^^ 21 — 1 — m, 22 >-^ 22 — ^ — n and using (19) one 

I Cone 3 hi{2i,22) 



W 



where 



hli2i,22) = Ul 



z\ — l~m Z2~l—n -p2 



U 



rp , — , . d2\ A d22 

Typei 2122(21 + 22) 

r(i + zi)r(i-zi) 



(62) 



r^(l+m-zi) 



T{2 + m-2i) 



X r' (1 + „-.,) ■- <i + -^> ■- <' - '-^> ^ (!,+ --■ + ^2) r (1 - .. - ..) ^ (^3, 



r(2 + n-Z2) 



r {2 + m + n — 2i — 22) 



We now fix the 2-vector / by the geometrical rule described when we presented the trans- 
formation law and choose the matrix A and 2-vector g as simple as possible. It is clear 
that the singular gamma functions of i?_i in this type of intersection are r(zi), r(z2) and 
T{1 + 2i + 22) and that the singular lines corresponding to the first and last of these three 
gamma functions intersect one side of the /-line, whereas the singular lines corresponding to 
r(z2) intersect the other side (see Figure [8|. Therefore, from the conventions of Figure Ig] one 
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\T 



finds / = [zi{zi + Z2), Z2) ■ It is easy to find that the simplest g and A which fulfill all the 
requirements of the transformation law are g = {zf, Z2)^ and 



Since det A 



1 one then finds 
Res 



A 



hi{zi,Z2) 



ZlZ2{zi + Z2) 



1 -Zi 

1 



Res 



hl{zi,Z2) 
zf Z2 



(64) 



(65) 



Therefore 




R- 


-l{ui,U2) 

00 00 


Cone 3 
Type 1 


= 


EE^^^o,o) 




m=0 n=0 






00 00 


,— m— 1„,— n— 1 



m=0 n=0 



Uo 



r(m + l)r(n + l) 



m + l)(n + l)r(2 + m + n) 



[In til - 2'i/'(l + m) + ■ip{2 + m) + V(2 + m + n)] 



t» 00 „.-ni-l„,-?i-l 



EE^ 



Un 



^r(m + l)r(n + l) 



m=0 n=0 



(m + l)(n + l)r(2 + m + n) 



1 



m + 1 



-f^m + Hi. 



■m+n 



+ Inui 



(66) 



Notice that if the choice of / is unique in the transformation law, this is of course not 
the case for A and g. But with other choices, one obviously finds the same result at the end. 
Indeed, instead of g^ = (zf , 22)"^ one could have chosen for instance g = (z^, ^2)^- In this case, 
one possible matrix of the transformation law would have been 



A 



Zl - Z2 
Z2 



Z1Z2 
1 - Z1Z2 - 



with det A = zi 



Z2- In this case one finds 



Res 



^^1(^1,^2) 

ZlZ2{zi + Z2) 



Res 



hi{zi,Z2) 
Zi Z2 



Res 



hi{zi,Z2) 
Z2 



Res 



hl{zi,Z2) 



(67) 



(68) 



By the Cauchy formula the second and third terms are zero and one obtains the same result 



as in eq. (66) 



Type 2: the singular points (21,^2) = (0, —1 — m) where m > is an integer. 



Performing the same change of variable for Z2 as in the case of Type 1 and no change of 
variable for zi one has 



UJ\ 



iCone 3 _ /?-2(^l,^2) 
^yP^ 2 - zlz2{zi+Z2) 



dzi A dz2 



(69) 



where 



h2{zi,Z2) = nf n^^-^~" T^ (1 - zi) T (1 + zi) T^{l + m- Z2) 

r (1 + Z2) r (1 - Z2) r (1 + zi + Z2) r (1 - zi - 22) 



T{2 + m-Z2) 



r (1 + m - zi - Z2) 



(70) 
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Looking at Figure pi (and Figure pi) one finds / = (^^i + Z2,ZiZ2)'^- Now, with g = (z^,^!) 
and 



A 



-1 



Z2{Z2 - Zi){zf + Z^) Zi 

whose determinant is det A = zf + zfz2 + Z2 — Z1Z2 — ^i^l' o^^^ finds 

h2{zi,Z2) 



(71) 



Res 

so that 



zfz2{zi + Z2) 



R-l{ui,U2)\ 



2 



^^^•^^(0,0) - \hf'\Q,Q) + ;Jr4''°^(0,0) - ;^4°''^(0, 0) , (72) 



3! 



3! 



Cone 3 
Type 2 



1 ^ tij'-'" 



E? 



6 ■'^— ' 1 + 771 
m=0 



hl^ + i/„ 

+ 



1 



1 + 771 
1 



In \- Hm- T— 

U2 1+777 



I8C2 + 



(l+m)2 



Olirn 



(l + 7n)3 
(we do not give here the equivalent expression in terms of polygamma functions). 



(73) 



Type 3: the singular points (zi, Z2) = (1 + ?77, —2 — m — n) where ?t7, n > are integers. 
Performing the change of variables zi i—?- zi + 1 + rri, Z2 '—^ -^2 — 2 — 7?7 — 77, one has 



UJ 



I Cone 3 
Type 3 



h3izi,Z2) 
zlz2{zi + Z2) 



dzi A dz2 



(74) 



where 



/73(Z1,Z2) = (-1)"*-^X 



m+1 zi+l+Ti Z2~2—m 



U. 



.„r2(i + zi)r2(i-zi; 



r(l + r77 + zi) 



2 r2(2 + 7n + zi) 

r (1 + 22) r (1 - Z2) r (1 + zi + za) r (i - zi - zs) 



X r^ (2 + 777 + 77 - Z2, , , , , 

r (3 + 771 + 77 — Z2) r (1 + 77 — Zl — Z2) 

Looking at Figure p^ one finds / = (zi + Z2, z^Z2) . Now, with g = {z\, z^) and 

A=( '" -1 

V 2:2(^:2 - Zl) 1 

whose determinant is det A = z^ + z| — Z1Z2, one finds 
hz{zi,Z2) 



(75) 



(76) 



Res 



ZiZ2(zi + Z2) 



i/7?^°)(0,0) + i/7f^)(0,0)-4^'^)(0,0), 



(77) 
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so that 



R-l(ui,U2 



I Cone 3 
iType 3 



Y, Y.{-ir^'u\ 



+m —2—m—n 

iir. 



m=0 n=0 



r(l +m)r2(2 + m + n) 
r2(2 + m)r(3 + m + n)T{l + n) 



Ul 



n2 



In — ^ + -ffm — 2//l+m + 2-fri+m+n — -f^' 



2+m+n 



-r<J',2 2 *" "T -"m+1 -"l+m+n "T 2 2+m+n 



(we do not give here the equivalent expression in terms of polygamma functions). 



(78) 



We numerically checked all these results as we did for the case of Cone 1. Notice that 
thanks to the symmetry of R-i it is straightforward to deduce the series representation 
corresponding to Cone 5 from the results of this section. 



3.1.5 Other cones, other regions 

As we saw in Figure HI there are two other cones associated to i?_i: Cone 2 and Cone 4 and, 
thanks again to the symmetry of (12) only one of them has to be considered. 



This is what we did by computing the corresponding series for Cone 4 and, once more, 
we checked numerically the expressions that we obtained by a comparison with the direct 



evaluation of the integral representation (12) for values of ui and U2 lying in the relevant 



convergence region 7^4. Since these results do not really bring important information, we do 
not list here the expressions of the series in order to shorten the paper. They are however as 
straightforward to obtain as the others. 

Nevertheless, it is perhaps instructive to briefly consider one particular type of contribu- 
tions, which belong to the overlapping regions of Cone 3 and Cone 4, to insist on one point 
concerning the derivation of the 2-vector /, in the transformation law. More precisely, we 
will have a look to the contribution of the singular points (zi, Z2) = {1 + m, —2 — m — n), 
where m,n > are integers, in Cone 4. One can indeed see that we already met this type of 
singularities in Cone 3 (Type 3). Therefore part of the calculation has already been done and 



one may directly take the result of eq. ( 74 ) and eq. ( 75 ) . 

The important point is that, after a look at Cone 4 in Figure 4 one finds / = (^2(^1 + 

Z2),z'f) which is not the same expression as the one of Cone 3, that we recall to be / = 
(zi + Z2,zfz2) (see the text after eq. (75)). This is of course due to the fact that the 
geometrical rule which allows to fix / does not give the same results in different cones even 
if the singular structure under study is the same in these cones. Now, with g = (z^ 



2 ^3)T ^^d 



A 



1 

Z2 - Zi Z2 



whose determinant is det A = zi — Z2, one finds 

h3{zi,Z2) 



Res 



Z(Z2{ZI + Z2) 



j4°'')(0,0)-4^'^)(0,0). 



(79) 



(80) 



where h^ is given in ( 75 ) . 
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Therefore, although in our notations the left hand side of (80) and (77) are formally the 



same, their right hand side are not, so that obviously the result given in eq. (78) do not 
coincide with the expression of R^-i(u-\,U2)\m ,^ . n %■ 

J- J- \ J-i '^'IType (1+m,— 2— m— n) 



3.2 The i?_i_z-^('Ui,-U2,M3) integral 

After our exposition of the method in detail for the simple case of the R-i integral considered 
in |8], we want to mention a few facts about the R-i-zj^{ui, U2, U3) integral of the same paper 
(second line of eq. (3.6) in [5]): 



R- 



1-2 



i{ui,U2,U3) 



c+ioo d+ioo 

dzi f dz2 



2'iri J 2'iri 

c—ioc d—ioo 



u\^u\-u-^^- 



X -r2 (-zi) r^ (zi + 1) r2 (-^2) r (z2 - ^1) r (zi + ^2 + 1) , (81) 
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where c 



-g and d 



1 

'4- 



The singular structure of the integrand of (81) is represented in Figure 9 



Rcc, 



-Re:i 




Figure 9: Singular structure of the integrand of R-\-z^{u\^U2iU-i) . The red dot is the point 
(c,d). 



Drawing the straight /-line and rotating around 7 = (c, d) it is easy to conclude that five 
cones are to be found, that we represented on Figure [TO] 

It is interesting to present (part of) the calculation of the i?_i_2-^(ni, n2, Us) integral for 
several reasons. 

At first sight it is a less simple integral than R-\{u\^U2) since one can see that it is not 
symmetric under an exchange of — and U2- Notice however that apart from its — term the 
integrand of (81) is symmetric under the replacement z\ — )• —z\ — 1 and this is reflected in a 
partial symmetry of the different cones (see Figure 10, where one can indeed remark that the 



23 




Rez2 



»- ReZ] 




Rezi 



Figure 10: The five different cones of R^i^zi{ui,U2,U3). Small black circles distinguish 
spurious singularities. 



symmetry between Cone 1 and Cone 2, as well as between Cone 4 and Cone 5, is broken only 
by the zi = line) and in fact also of their corresponding regions of convergence. 

The calculations of the different cone's contributions are not harder than the computations 
of previous sections and one may follow the same systematic evaluation method: we do not 
need to perform any manipulation at the integral level (by inserting Euler representations) as 
the authors of |71 18| did (see Appendix C.2 of [8j). 

Moreover, one will see that the regions of convergence of the different series representations 



of (81) are slightly more complicated than those of R-i{ui,U2) (see Figure 11). This means 



that the region of convergence of the series representation of the two-loop hexagon Wilson 
loop (see [311]) is probably tricky to find (without the help of analytic continuations). 

At last, the R^i^ziiui,U2,U3) integral presents the feature that (as we already said in 
section 3.1 ), for two of its cones, some subsets of singularities are spurious. This fact, that we 
did not meet in the case of R-i{ui,U2), should however not be thought of a rare occurence 
in general and we think that it is important to show explicitly on this simple example how 
spurious singularities may appear. 



3.2.1 Regions of convergence 

Let us now give the region of convergence of each cone (to obtain them we follow the same 



procedure as in section 3.1.2). 



Conel The region of convergence is given by TZi = {\u2\ < 1 and 



+ 2 



(yellow region in Figure 11 ). 



UlU2 
"3 



<1} 



Cone2 The region of convergence is 7^2 = {|''^2| < 1 and 
region) . 



+ 2 



"1 



< 1} (green 
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Figure 11: Regions of convergence of the series representations of R-i-zi{ui, U2, U3). 



ConeS The region of convergence is TZ-^ 
- } (cyan region). 



"1 



M2M3 



1 + A 1 



< 1 + Wl 



and 



M3 



«2«1 



< 



Cone4 The region of convergence is TZ^ 
< 1} (orange region). 



UIU2 

U3 



U2U3 



< 1 + Jl 



and 



+ 



ConeS The region of convergence is TZ^ = { 



U3U2 



< 1} (magenta region). 



< 1 + Jl 



and 



+ 



The (partial) symmetry mentioned above between the cones is clearly seen here as we 
perform the change i^ — ^ ^ and corresponds indeed to the change zi — )• — zi — 1. It is 
however completely manifest at this level since the — term of the integrand has no influence 
on the different regions of convergence. 



As for the case of R-i{ui,U2) one can see in Figure 11 that an infinite white band is not 



reachable by any of the series representations that our method allows to derive. 

Contrary to the case of i?_i(ni,U2) the regions of convergence T^i, ..., 7^5 cannot be 
obtained from the results of section |B.5[ since the corresponding series are not Kampe de 
Feriet double series. 

We give now the results of the calculation of the series associated to Cone 1. 
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3.2.2 Cone 1: region of convergence IZi = {|n2| < 1 and 



+ 2 



UlU2 
"3 



<1} 



As can be seen from Figure 10 tliere are 5 different types of intersections in this cone. 



Type 1: the singular points (21,22) = {1 + m,2 + m + n) where m,n >0 are integers. 
This type of singularities does not require to use the transformation law and one finds 

n / N I Cone 1 

'+" .2+n+n^^(l + n)^(4 + n + 2m) 



EE 

m=0 n=0 



Uo 



(m + l)r^(r7T, + n + 3) 



In U2 - 2'ip{m + n + 3) + ^p{l + n) + ■ip{4: + n + 2m) 



In — - ij(m + 2) + V'fm. + 1) - ^(1 + n) + V'(4 + n + 2m) 
U3 



^^^\l + n)+ ^(^)(4 + n + 2m) 



(82) 



Type 2: the singular points (zi, Z2) = {1 + m, 0) where ?ti > is an integer. 
In this case the transformation law is necessary and one finds 

J-, I \ I Cone 1 

00 



m=0 



Type 2 
Ui \' ' '^ 1 
U3 



m, + 1 



til 

In h -0(1 + m) — 0^(2 + m) 



("in— +-0(l + m) - V'(2 + m)') ( - In — + --0(1 + m) + -lnn2 - ipil) + --0(2 + m' 

\ U3 J yp „. c n a ' 



6 M3 6 



+ ^(^^(1) + ^V'^'^(l + rn) +^0«(2 + m) 



+ 



lnii2 - 2V'(1) + 2V'(2 + m) 



3V'^^^(1) + ^-^^^\\ + m) - -0W(2 + m) 



+ ^^(2^(1 + m) + ^V'^^^(2 + m) I. 

O D I 



(83) 



Type 3: the singular points (21,22) = (1 + m + n,l + n) where m,n >0 are integers. 
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This type of singularities does not impose to use the transformation law and one finds 

I Cone 1 



R-l-ziiui,U2,U3) 



Types 



EE(-i) 



m=On=0 



IH-m+n 
U3 



r(3 + m + 2n) 



(1 + m + n)r(l + m)r2(2 + n) 



1 
3! 



Ul 



In — ^ — V'(2 + m + n) + ^(1 + m + n) — -0(1 + tti) + ijj(3 + 2n + m) 
^3 



X 



In 7/;(2 + m + n) + ipil + m + n) + 3 In n2 

^^3 



6^(2 + n) + 2V'(1 + m) + 4V'(3 + m + 2n) 



+ 3 



Ul 



+ 6 



In — ^ - ^^(2 + m + n) + ^(1 + m + n) + In n2 - 2V'(2 + n) + 2V'(3 + m + 2n) 
^3 

6V'^^^(1) - V'*-^-'(2 + m + n) + tp'^^^l + m + n) - t/j'-^Xl + m) + tp'^'^^S + m + 2n) 
In -0(2 + 771 + n) + ^(1 + 771 + n) — -(/'(l + m) + 7/>(3 + 2n + m) 

U3 
X 



-2V'^^Hl) + V'^^^l + m) + ^(i)(3 + m + 2ra) 
V^(^)(2 + m + n)+V*^^^(l + m + n) + 2V'^^^(l + m) + 4V'^^^(3 + m + 2n) \. (84) 



Type 4: the singular points (21,^2) = (0, 1 + m) where m > is an integer. 

This type of singularities does not require to use the transformation law and one finds 



Ti I \\ Cone 1 



1 



E< 

m=0 



1+m 



1 



?Tt + 1 



lnu2 - 7/;(2 + m) + 7/^(1 + m)\ ^^/'(^^(l) + V^^^l + 



m) 



+V^(^)(2 + m) + ( In — - V(l + rra) + V(2 + ra) 

\ Us 



xf In— -ipil + m) + 'ip{2 + m)] + i{j^'^\l + m) + V^^^(2 + m) 

Type 5: the singular point (^1,^2) = (0,0). 

For this last point the transformation law is necessary and one finds 



R-l-ziiui,U2,U3)\ 



Cone 1 
Type 5 



V'(')(l)-[3V'«(l)+^ln2'^ 
L 6 U3 



lnti2 



2 ^ 24 us 



ln"i ln"i 



U3 J Us 



(85) 



2 U^'^^ (1 + m) - V^^^ (2 + m)) (86) 



(87) 



^[V.«(l)]^-l^(3)(i). (88) 
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3.2.3 Other cones 

As for the caseoi R^i{ui,U2), we have computed the series representations oi R-i-^iiui, U2,U3) 
coining from all the cones of Figure [TO] and checked that the results were in agreement with 



direct numerical evaluations of the integral representation (81) for some values of ui, U2 and 
^3 in the corresponding different regions of convergence. 

The contributions of all cones were easy to compute but we will not write the corresponding 
results to shorten the paper. 

One point has however to be underlined: in each of the cones 4 and 5 there is an infinite 
subset of singularities which do not have to be included in the calculation of their respective 
contributions. For Cone 4 these singularities are located at [zi,Z2) = {2+n, — | —m — n) and 
for cone 5 they are at (zi, Z2) = (— | — m, — | — m — n), where as usual m,n > are integers. 
We recall that the reason why one has to discard these singularities is that in each case all 
singular lines which cross each other at these singular points are also crossing the same side 
of the /-line corresponding to the cone under consideration. Therefore they do not fulfill the 



condition presented in section 3.1 for being considered as relevant singularities. Notice that 
these singularities which are spurious for Cone 4 and Cone 5 are however completely relevant 
for the calculation of the series associated to Cone 3. 

3.3 The scalar box integral with one external mass 

As a last twofold example, we consider the scalar box integral with one external mass. This 
divergent integral, computed in dimensional regularisation, will give an example of how the e 
singularity may be simply solved graphically with our approach. 

This integral is also a nice example since, due to its simplicity, its double series represen- 
tations may be resummed in terms of simple functions as well as combinations of the Gauss 
hypergeometric function 2F1. By well-known analytic continuation formulas, this allows to 
reach in a simple way the regions in the parameters space which are not directly accessible 
by the series representations. 

Modulo an overall factor, the box integral with one external mass in D = A+2e dimensions 
has been expressed in |15) as 



c—ioo d—ioo 

xT{2-e + zi + Z2)T{e - 1 - zi - Z2)T{1 + Z2)r(e - 1 - Z2)T{1 + zi + Z2) , (89) 
where the conditions c<0, d — e<— 1 and c + d > —1 follow from the integration over 



Feynman parameters and fix the fundamental polyhedron (the red triangle in Figure 12). 



The singular structure of the integrand of (89) is represented in Figure 12 



It is obvious from Figure 12 that the e dependence of the integral implies a splitting of 
some of the singular lines that would be identical if e were zero and, therefore, that this will 
allow to obtain the results in a simpler form since the multiplicities of corresponding the poles 
are decreased. It is however clear that when e — t- the fundamental triangle collapses, giving 
birth to a pinch singularity located at (zi, Z2) = (0, —1) which of course reflects the divergence 
of the integral. 

Drawing the straight /-line and performing the rotation around 7 = (c, d), it is easy to 
conclude that only five cones may be found, that are represented on Figure [13] Cone 2 
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Figure 12: Singular structure of the integrand of the box integral (the right figure is a zoom 
of the left one, in order to better show the fundamental triangle). The black dot is the point 
7 = (c,d). 



corresponds in fact to the calculation performed recently in |15] . which allows a cross-check 
of the calculation. Notice also that for each of the three integrals considered until here in 
this paper there were five differents cones. However, this is of course not always the case in 
general. 

The series representations of the box integral associated to all the cones but Cone 2 do have 
polygamma functions in their general term, even if e is kept at a non-zero value, since they 
contain subsets of singular points coming from the intersection of three singular lines. The 
simplest cone is therefore Cone 2 whose singular points are formed only by the intersection of 
two different singular lines. In fact, its corresponding double series may be resummed in terms 
of the Gauss hypergeometric function [15j and this allows, since one knows the e expansion 
of this function |T^, to obtain the e expansion of the box integral from an exact result. We 
will compare the first few terms of this expansion with the result that may be obtained with 
our graphical approach for solving the e singularities. 

Before this we discuss the convergence properties of the series representations of the box 
integral, and we give the expressions of the series associated to Cone 1 and Cone 2 (the series 
corresponding to the three other cones have also been obtained easily). 



3.3.1 Regions of convergence 

From Horn's theorem, one gets the following regions of convergence (summarized in Figure 
141): 



Cone 1 The region of convergence is given by TZi = { 



(magenta region in Figure 14). 



+ \j\ < 1 and 



< 1} 



Cone 2 The region of convergence is 7^2 = {X + XX-'^J' (orange region). 
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Figure 13: The five different cones of tlie box integral. Tlie daslied lines are those which 
have an e dependence and the small black filled circles distinguish spurious singularities as 
described in section [3.1.21 



Cone 3 The region of convergence is TZ^ ~ { ^ + s^ ^ ^^^ ~^ < 1} (yellow region). 



Cone 4 The region of convergence is 7^4 = { ^ + If I ^ ^ ^^'^ 1^1 ^ ^} (cyan region) 



Cone 5 The region of convergence is TZ^ = { 



+ - < 1 and 



< 1} (green region). 



However, as we will see explicitly in sectiom ^'^||3.3.2 and 3.3.3 for the cases of Cone 1 and 
Cone 2, the double series representations of the box integral may be reduced to combinations 
of the Gauss hypergeometric function 2F1 and of some trivial functions. Therefore, well- 
known analytic continuation formulas allow to cover the whole ( I -^ I , 1^1) plane but a few 
exceptional singular lines. 



All results given in these sections have been checked by a comparison with the direct numerical computa- 



tion of ( 89 1 for particular values of the expansion parameters and e. 
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tls2 




12 3 4 

Figure 14: Regions of convergence of the series representations of the box integral obtained 
from Horn's theorem. 



To be more precise, we give the following simple concrete example, taken from ||9j, which 



is also met here in ( 92 ) and ( 99 ) . 



Let us consider the series 



oo oo 






r?i=0 m=0 



m\ n\ 



(90) 



The region of convergence of this series is given in ^j as 7^ = |ti| + |t2| < 1, which is also 
what one may conclude from Horn's theorem. 



However one may rewrite ( 90 ) as 



^(ti,t2) = F2(a,/3,/3',/3,/3';-ti,-t2) 



(91) 



(l + tl+t2)"' 

where F2 is one of the Appell hypergeometric series. It is clear that the domain of definition 



of the right hand side of (91) is not restricted to IZ, and allows an analytic continuation on 



the whole (|ti|, \t2\) plane except on the line 1 + ti + ^2 = 0. 

3.3.2 Results for Cone 1 

Notice that in this cone there are two different subsets of spurious singularities located at 
(2^1, Z2) = {1 + k + n, —1 — k) and (zi, Z2) = {e + k,n), where k,n > are integers. 
Let us now give the results corresponding to the relevant singularities. 
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Type 1: the singular points (zi, Z2) = {k,n — 1 + e) where k,n > are integers. 



I{s,t, m ) 



Cone 1 



Typel (47r)2+': yfj,"^ 



tV i-t)-^ ft 



1-e 



vr 



r(2e) \s J sin(7re) 
00 00 ^ 2 \ '^ 



k=On=0 
t 



m 
T 



s\" r{n + k + e) 
t) r(l + A:)r(l + n) 



In - + ^(1 + n) + V'(l - e) - V'(e) - V'l^ + A: + e) 
s 

r(2e) \s) sin(^e) ^^^ V t ^ t 

ln-+V'(l-e)-V'(6) + -f-l-- + — y 2Fif/'' 



(47r)2+^ V/^^ 



1 + 



, , m \ / s 



In -1 + 



m"^ t 



s s 
(92) 



Type 2: the singular points (21,^2) = {k,n) where k,n>0 are integers. 



I{s, t,m 

00 00 

<EE 

k=0 n=0 



I Cone 1 



tV(-i)- 



\Type2 (47r)2+^ V /^ V r(2e) 

m^y / sY r{l + k + n)r{2 + k + n - e)r{-l - n + e)r{-l - k - n + e) 



t 



r(l + A;) 



(93) 



Using ( 19 ) this double series may be expressed as a special case of the Appell hypergeometric 
function F2 which in fact may be reduced to the Gauss hypergeometric function 2F1 



I{s, t,m'^) 



Cone 1 



-t 



i-ty 



vr 



Type2 (47r)2+^V/^V T (2e) T (2 - e) Vsin(7re) 

t 



t-m2 2 H2-e 



m^ 



1 . (94) 



Type 3: the singular points (zi, 22) = {e + k + n, —1 — n) where k,n > are integers. 
Conei i /_tY {-t)-^ fm^Y t 



I{s, t,m 



Types (47r)2+^ \fj.'2j T{2e) \ t 



m 
t 



2\ '^ 



m 



2\ " 



T{-k-n-e)T{k + e)T{n + e). (95) 



xEE 

fc=0 n=0 

Once again this series may be reduced to the 2F1 function with the help of the generalized 
reflexion formula. It is a special case of the Appell function F3 which can be expressed as a 
special case of the Appell function Fi, and the latter simplifies to the 2F1 function as follows: 



I{s, t,m ) 



Cone 1 



Type 3 



(4^)2+^ V A*^ 



-m 



i-t)-' vr r(e) t 
r (2e) sin(7re) e s 



2-ri 



1 + e 



l-(l-"^)(l-^ 



(96) 
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Type 4: the singular points (zi, Z2) = {1 + k + n, —2 — k + e) where k,n > are integers. 



I{s,t,m 



Cone 1 



tV {-t)-^m^ ft 



Type A (47r)2+^ V /^ V r(2e) t \sj sin(7re 



2-e 



TT 



fc=0 n=0 



2\ 'i: 



m 



2\ " 



r(fc + l)r(n + e) 
r{k + n + 2) 



(97) 



In this case, one first finds that (97) may be written as a special case of the Appell function 



F3 which can be expressed as a special case of the Appell function Fi, and the latter simplifies 
to a sum of 2F1 functions as follows: 

Cone 1 ■ / ^\ e / ^\-2 _ /a\ 2-e 

Type 4 



I{s,t,'m ) 



(4^)2+^ \fl^ 



ty {-t)-^S ft 



TT Tie) t - m^ 



2^1 



1,1 -e 
2-e 



r(2e) t \s J sin(7re) e — 1 m? — s — t \t — m? 



t 



1 



s + t — m^ 



+ 



t 



t — m? 



2-ri 



1,1-e 
2-e 



t 



s + t — m^ 



(98) 



3.3.3 Results for Cone 2 

As we said above, the calculation corresponding to this cone has been considered recently 
in |15) . We give anyway the final results since it will be useful in the following. Indeed, 
this cone allows for the simplest exposition of the way we treat the e singularities with our 
computational approach. 

The double series corresponding to the following four different intersections types are 
special cases of the Appell F2 function which can each time be simplified. 



Type 1: the singular points (zi, Z2) = (— n — e — k,n — 1 + e) where k,n > are integers. 

Cone 2 -' / + \ <= /' +^-2 -i . _ . , / _ \ 2 



I{s,t,m 



t V(-t)-2f 



vr 



Typei (47r)2+^ ^ fi"^ J T{2e) s \m?J Vsin(7re) 

T{k + n + e) 



EE 

A:=0 n=0 
I I t 



t 



m^ 



(4vr: 



i2+e 



/"^ 



m^ 



1 1 



r(l+n)r(l + A;) 



r(2e) st V sin(7re 



TT 



r(e) 



TTi^ — s — t 



(99) 



Type 2: the singular points (zi, Z2) = ( — 1 — n — k,n) where k,n > are integers. 



I{s,t,m ) 



Cone 2 



t Yi-t)-^ t 



Type2 (47r)2+'s ^ A*^ / r(2e) Tm? \sin(7re 



vr 



EE 

fc=0 n=0 

i 

(47r)2+^ 
1 



?n,^ 



m 



r(fc + n+l) 



s \" 

~^! T{l + k)T{2 + n-e) 



^^ 



1 



TT 



r(2e) Vsin(7re)y r(2 - e) 



1 



"" t(m2 - t) 2^H 2 - e 



m^ — t 



(100) 
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Type 3: the singular points (zi, Z2) = (—2 — k — n + e,n) where k,n > are integers. 



I{s,t,m 



Cone 2 



tVi-t) 



Type3 (47r)2+^ V MV r(2e) \m'^ 



-2 



2-e 



vr 



sin(7re) 



EE 

fc=0 n=0 

i 



t y / s \" r(A: + n + 2-e) 
^y V^J r{k + 2-e)r{n + 2-e) 

2 



2 \ f 1 
m^\ 1 



(4^)2+^ V /"V r(2e) Vsin(7re)y r(2 - e) 



vr 



1 



1 



m 



2 — f m^ 



^ ^ ' 2-e 



St 



{m? — t){m'^ ~ *)/ 



(101) 



Type 4: the singular points (zi, 2:2) = ( — 1 — n — /c,n— 1 + e) where k,n >0 are integers. 



I{s,t,m 



Cone 2 



Type A 



fJ- 



t_y i-t)-^ t ft 

~2 



T{2e) m^ \ s 



l~e 



vr 



sin(7re) 



s \" 



r(/c + n+l) 



(4vr)2+e 

00 00 

EE 

fc=0 n=0 

i 
(47r)2+^ V AiV r(2e) Vsin(7re)7 r(2 - e 

"^' / 1,1 

TTT,^ — s \ z — e m^ 



t 
m? ) \m?J r(l+n)r(2 + A;-e) 

2 



1 



TT 



1 



(102) 



3.3.4 Cone 2: e expansion 

Until here our results have an exact dependence in the e dimensional regularisation parameter. 



We will now show how to solve the e singularities in a simple way from Figure 12 , for the case 
of Cone 2 (obviously the same method may be used for the other cones). 

Let us first recall that it is easy to compute the first few terms in the e expansion of the 
exact results obtained in section 3.3.3[ by using the e expansion of the Gauss hypergeometric 
function fl61 



2F1 



l,l±e 

2±e 



lie 

e->0 Z 



[- ln(l -z)^e Lis {z) + <? Lis {z) + 0{i')\ , 



(103) 



and one finds that the divergent part of the result associated to Cone 2 is of the form 
where 



A_2 



+- 



A. 



M^st 



and 



A. 



8-K^st 



jE — ln(47r) + In 



St 



m'^fj? 



(104) 



(105) 



In fact, without using the exact result as a starting point, it is very simple to get these terms 
(z.e to solve the e singularities), just by having a look at the picture representing the singular 
structure of the integral, and the one where Cone 2 is shown. Indeed since the fundamental 



triangle collapses when e goes to zero (see Figure 12), to avoid the pinch singularity it is 



necessary to "move" the 7 point outside the fundamental triangle in order to put it in a safe 
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region. There are many ways to do this and, as an example, we show two different possibihties 



in Figure 15, where one can see that the new 7 points (71 and 72) will not be pinched if e is 
taken equal to zero. 
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Figure 15: Solving the e singularities. The black dot is the point 71 = (c\^d\) = ( — 1,-2) 
the red one is 72 = (02,^2) = {-\,-\)- 



2^' 



Now, if one looks for the cones associated to the integral I{s, t, m?) with 7 replaced by 71 
(resp. 72), one finds for the relevant conc^^ the region C2 (resp C2) (see Figure 15). It is then 
obvious that 

I{s,t,m^)\?;-^'^ = I{s,t,m')\^^^ + Res.i (106) 



or 



/(s,t,m2)|C°-2^/(s,t,m2)|g+iies.2, 



(107) 



where 



Res.i = ^ [-Res.|(_,_fc _i+,) + i?es.|(_i_fc _i+,) + Res.\(_i_k,k) + -Res.|(_i_e_A;,fc+,)] (108) 

fe=0 
and 

CX) 

Res.2 = ^ [-Res.|(_,_fc _i+,) + Res.y_i_k-i+e) - Res.\(^_i_k^k+e) - -Res.|(_i_,_fc^fc)] (109) 

fe=0 

and that the -j and - terms will of course be entirely given by Res.i or by Res.2 (which we 
checked explicitly in both cases). 

4 Mellin-Barnes integrals of higher dimension 

In this section we give a foretaste of the extension of the method to Mellin-Barnes integrals 
of higher dimension. We only consider here a simple threefold integral, as a toy model, and 
one of its n-dimensional extension, since the general procedure is presently under study a will 
be treated in a forthcoming publication. 



^^By relevant cone we mean the cone which has the largest overlap with Cone 2. 
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4.1 A toy integral 

Let us consider the integral 



I{ui,U2,U3) 



c+ioo d+ioo e+ioo 

ds f dt f du 



2m 



27ri I 27ri 



c—ioc d—ioo 



X ui u\ -u^ r^ {-s) r^ {-t) r^ (-u) r^ (i + s + 1 + -a) , (no) 



where c 



d 



-g and e 



10- 



As before, to obtain the triple series representations of (110) we have to find the different 



cones associated to this integral, in order to know the different convergence regions as well as 
over which residues the corresponding triple sums will be performed in each case. However, 



the singular structure of the integrand in (110) can not be represented in a picture as readable 



as, for instance. Figure [9] and, moreover, we do not have a simple geometrical procedure to 
obtain the cones as the one used in the twofold case. Nevertheless it is possible to obtain 
them from an exhaustive search of the possible simultaneous constraints coming from the 
different triplets of gamma functions, in the same spirit as what was done in the beginning of 
section 3.1.1 Indeed, only the following four different regions of the (Re s,Re i, Re u) three 



dimensional space have to be considered: 



Cone 1 
Cone 2 
Cone 3 
Cone 4 



Re s > 0, Re t > and Re n > 0. 
Re s > 0, Re t > and Re (1 + s + t + u) < 0. 
Re s > 0, Re u > and Re (1 + s + t + u) < 0. 
Re t > 0, Re n > and Re (1 + s + t + u) < 0. 



As can be seen on Figure 16 the cones are in fact infinite pyramids, the fundamental 



polyhedron of the integral, given by the constraint Re s < 0, Re t < 0, Re n < and 
Re {1 + s + t + u) > 0, being a tetrahedron. 

It is easy to see that the singularities in Cone 1 are only of Cauchy type, therefore the 
corresponding residues will be trivial to compute. The three other cones, however, will neces- 
sitate a three dimensional extension of the transformation law. In fact, due to the symmetry 
of the integral, one will see that in principle only one cone has to be computed, the three 
others being deduced from the latter by simple changes of variables. This will give us the 
opportunity to check that our implementation of the transformation law at the 3-dimensional 
level is correct. 



4.1.1 Cone 1 

Let us begin by the calculation of the triple series corresponding to Cone 1. In this simple 
example there is only one type of intersection in the whole cone: the points (s, t, u) = (m, n,p) 
where m,n,p > are integers. 



has 



Performing the change of variables s >-^ s + m,t >-^ t + n,u i-^ u + p and using ( 19 ) one 

h{s,t,u) 



UJ 



s^t'^v? 



ds /\dt /\ du 



(111) 
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Figure 16: The fundamental tetrahedron and different cones of I{ui, U2, U3). 
where 

,. , .^ s+m t+n u+p T^ [1 + S)T\1- S) T\l + t)V^l-t) 

n{s,t,u)-u^ u^ ug p2(^ + ^ + i) r2(t + n + l) 

r2(i + ^)r2(i-n)^2 



r2(« + p + i) 



r^(s + t + u + m + n + p+ 1) . (112) 



Using now the three dimensional extension of Cauchy Formula (22) one has 

5^/i(s,i,n) 



I{u\,U2,Uz) 



I Cone 1 



00 CXD 00 

EEE 

m=0 n=0 p=0 



ds dt du 



(113) 



(0,0,0) 



the overall sign coming from the ( — 1)^ due to the closing of the integration contours. 
It is then straightforward to show that 



I{ui,U2,Uz)\ 



Cone 1 

00 00 00 



EEE 

m,=0 n=0 p=0 



n^n^n^ 



r2(l + m + n + p) 

r2(i + m)r2(i + n)r2(i + p) 



X i f{ui,m + n + p,'m)f{u2, m + n + p,n)f{u3, m + n + p,p) 

+ 2ip^ '{1 +m + n + p) [/(ni,m + n + p,m) + f{u2,m + n + p,n) + f{uz^m + n +p,p)] 

+ 2^(2)(l + m + n + p)|. (114) 

where 

f{x,y,z) = \nx + 2ij{l + y)-2i){l + z). (115) 
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From the results in sections B.4 and B.3 this triple sum converges in the region 
U2\ + y^lual < 1 (notice that this region of convergence may also be obtained from the 



statements of section |B.5 since the series with general term /i(0, 0, 0) is a triple Kampe de 
Feriet series). 

For particular numerical values of the parameters ui , U2 and U3 in the region of convergence 
we checked that one recovers the value obtained from a direct numerical evaluation of the 



integral (110), with the same values of the parameters, to a very good accuracy. 



4.1.2 Cone 3 



Once more, in the simple example of (110) there is only one type of intersection in the whole 
Cone 3: the points (s, t, u) = (m, —1 — m — n — p,n) where m,n,p > are integers. 

Performing the change of variables s >-^ s + m,t >-^ t — 1 —m — n—p,u >-^ u + n and using 



(19) one has 



CO 



h{s,t,u) 



s^v?{s + t + uY 



ds A dt A du 



where 



h{s,t,u) = u 



s+m ^t-1-m-n-p ^ u+n 
1 



U. 



Un 



r2(i + s)r2(i- 

r2(s + m + l) 



T"^ (l + m + n + p-t) 



r2(n + n+l) 



r {s + t + u + m + n + p + 1) 



(116) 



We now have to apply the transformation law, therefore we have to fix the 3-vector / and 
choose a 3 X 3 matrix A and a 3-vector g as simple as possible, in order to use the natural 



three dimensional extension of (23). In our present case of study the 3-vector / may be 



obtained without a general procedure as the one presented for the twofold case since there are 
only three different sources of divergent behaviour: t, -^ and . , , , — yr. Therefore, one finds 

/ = (s^jU^, {s + t + u)^) . Although not completely straightforward, it is also not too hard 
to find a relevant matrix A and its corresponding vector g. For instance g = (s^, t^, u^) and 
the matrix elements 



^11 = A32 = 1 

A12 = Ai3 = A31 = A33 = 

A21 = 2st{u + A) + 5^(3 + t) + 2*2(3 + s + u)-6u'^ + t^, 

A22 

A23 



-2tu{s -4) + 2*2(3 -s-u)- u^{t - 3) - t^, 
-2s{t - 3u) - 2tu - 5^(3 + t) + f + u^{t - 3). 



(117) 



fullfill all the requirements of the transformation law. 

Since det A = 3s^ + 2st + s^t — t^ — 6su + 2tu + Su^ — tv?, one then finds, keeping only 
the non-zero contributions 



Res 



h{s,t,u) 



s'^v?{s + t + uY" 



-Res 



6h{s,t,u) 



+ Res 
-FRes 



2h{s,t, u) 
2h{s,t, u) 



V 



u 



Res 



h{s,t,u) 



(118) 
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so that 

I{ui,U2,U3)\ 



oo oo oo 

Cone 3 / > / ^ / > 

m=0 n=Op=0 



d^h{s,t,u) 



di^ 



+ 



d^h{s,t, 



u 



+ 



(0,0,0) 

d^h{s,t,u) 



dt"^ d 



(0,0,0) 

d^h{s,t,u) 



(0,0,0) 



ds dt du 



ds dt^ 
Using the same way to write the derivatives than in section 2, one has 

I{ui,U2,U3) 



:il9) 



(0,0,0) 



I Cone 3 

oo oo oo 



EEE 



m^ — 1— m— n— p n 



Ui U2 



U3 



r2(l + m + n+p) 



m=0 n=0 p=0 



r2(i + m)r2(i + n)r2(i + p) 

m + n + p,m ] f [ — ,m + n + p,n] f I -^,7n + n + p,p 



+ 2'tp^'^\l + m + n + p) 



\U2 



\U2 



f { — ,m + n + p, m) +/| — ,m + n + p,n 

\U2 ) \U2 



1 



+/ ( — ,m + n + p,p 

\U2 



+ 2'ip'^^'> {1 + m + n + p) } , (120) 



ui 



:, U2 



— and U3 — > ^. 

M2 M2 



where f{x,y,z) is given in (115) and which is nothing but (114) multiplied by — , with 



The triple sum ( 120 ) converges in the region 



+ 



< 1. 



Once more, we mention that all results have been checked numerically. 



Moreover, we see that in the simple case of (110), the triple series representation corre- 



sponding to Cone 3 could have been derived from the one obtained from Cone 1. This is of 
course due to the fact that under the change of variable s ^ —1 — s — t — u we have 



I{ui,U2,U3) 



— T ( !^ J_ '^ 

U2 V ^2 ' W2 ' U2 



:i21) 



and it provides a check that our implementation of the 3-dimensional transformation law is 
correct. 

Thanks to this nice symmetry of the integral, which may also be observed at the n- 



dimensional level in (122), it is then possible de get the n-dimensional version of (120) directly 



from eq. (129), by mutiplying it by — and by doing a replacement of the arguments similar 
as the one we did above. 



4.2 Extension to n dimensions 

Notice that, if not totally obvious, it is however possible to get the n-dimensional extension 



of (114) corresponding to the integral 



Cl+lOO 



I{ui,...,Un) 



dsi 
27r~i 



Cl—lOO 



c„+ioo 



dSn 

2m 



X Ul' • • • <" rP(-Sl) • • • TPi-Sn) TP{l + Sl + --- + Sn) 



fl22) 
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where the Cj [j = 1, ...,n) belong to the n-dimensional fundamental hyperpolyhedron (which 



reduces to the red tetrahedron of Figure 16 when n = 3) and p is an arbitrary positive integer. 
Indeed, the residues to compute in this case, in the n-dimensional equivalent of Cone 1 
defined by Re si > 0, ..., Re Sn > 0, are those of the singular points (si, ..., s„) = (Mi, ...,Mn) 
where Mi, ...,Mn > are integers. Performing the change of variables Sj i— )• Sj + Mj (j = 
l,...,n) and using (19) one has 



where 

h{si,...,Sr. 



u 



,«l+A^l . . . „«n+Af„ 



h{si,...,Sn) 

P P 

^l ' ' ' ^n 



dsi A • • • A dsr. 



(123) 



n 



•n„ 



rp(i + si)rp(i-si) 



X • • • X 



TP (1 + Sn) TP (1 - Sn) 



rP(si+Mi + l) ■■ TP{Sn + Mn + l) 

X r^ (1 + Ml + • • • + M„ + Si + • • • + s„) . (124) 



The n-dimensional Cauchy formula then allows to write 

1 9"(P-i)/i(si,...,s„) 



^(«l'->^")|conel = (-irE---E 

Mi=0 M„=0 



(p-l)!« ds\-^ ■ ■ ■ dsl-^ 



(125) 



(o,...,o) 



To compute the multiple partial derivative, one needs the following relation [Tj 



QkiA \-k„ ^=1 '^n /^ 

^^^7t::^^m = E---Eu 



nj dx\' ••• dx'n" dxl'-'' ■■■ dx^rl 



^Kfl ^'. 



n ^n 



(126) 



as well as the Faa di Bruno formula 



d'' 



dx 



k •' ^"' / > W'm-l . 



k\ 



■ kl'^'^mil ■ ■ ■ rrikl 



f 



{rriiA hnife) 






dPy 
dxJ 



(127) 



where the sum runs over all fc-tuples (rrii, . . . ,mk) of non-negative integers satisfying the 
constraint mi -|- 2m,2 + 3m,3 -!-••• + knif^ = k, and the Leibniz rule 



d'^ \ _ V^ f^\ ^^"^ '^^ ^^ 
^(««) = 2^[J^e- -^g—r 



i=Q 



(128) 



With the help of these relations, the calculation of (125), if a bit heavy, is however easy to 
perform and one finally finds 



I{ui,...,Un)\ 



(-1) 



CXD p— 1 p— 1 



Cone 1 



{P 



1)1" Z^ ' ' ' 2-^ 2-^ ' ' ' 2-^ \ ^-1 

'' Afi=0 A/„=OA:i=0 A;„=0 ^ 



p — 1 



p — 1 
h 



n 



i=i 



Af, 



'^ ' ^A;,- - / 



r(M, + i)2Ei/T"'^^gi^ jpi!--.. 



n(-2V'(*-i)(M,- + l))"" 



kj~l~q 

E 

r=0 



Kj i 



.<^ llmi . . . 



/I^^mi! • • • TUql 



ri 



r!'"''mi! • • • m^! 
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L 

-l-q 

n 






y 



{kj - 1 



{k 



2^ 



(i-i)/ 



1) 



«=i 



,j - I - q - ry.""'^ ' " '■mi!---mfc^_/_g_^! 



Al 



1 + E^^^.- Elm^rTTT 



i=i 



Al^-^mil-'-mA! 



mf 






j=l 



(129) 



where A = n(p — 1) — Yll=i ^j ^^'^ where the sums without indices are understood as for 
(1271). 



Although the expression (129) looks a bit comphcated, it is straightforward to extract 



from it simple formulas like (114) for particular values of p and n. 



More generally, the results corresponding to the two other cones of Figure 16 as well as 



those of the n — 2 other cones of the n-dimensional integral (122), may be straightforwardly 



obtained from (114) and its n-dimensional version (129) by simply exchanging the arguments 

— ,..., Un —)• ^- Of course one also has to multiply by an overall 



as follows: ui 
—- factor. 






5 Conclusions 

In this paper, we have shown with a general method how to derive, from a given twofold 
MB integral, several convergent series representations as well as their associated regions of 
convergence, the latter being obtainable before a full computation of the general term of the 
series. MB integrals of higher dimension have also been considered. 

In the twofold case, the method may be based on a very simple graphical approach: once 
the singular structure of the integrand has been drawn, it is easy to perform the extraction of 
different sets of poles whose residues give the mathematical expressions of the different series. 
The extension of the method to MB integrals of higher dimension is possible in principle 
although, obviously, one cannot keep a graphical view, at least when the dimension is greater 
than 3. 

In each of the examples treated for illustration, we began by looking for a parametrization 
of the relevant sets of singular points and deduced straightforwardly the corresponding regions 
of convergence of the associated series. Then we computed the residues using simple results 
of multidimensional complex analysis. 

The different series are analytic continuations of one another and an interesting point of 
the method is that, as we have shown, one does not need to perform the resummation of 
the series corresponding to a given region of convergence and to use analytic continuation 
properties of the result, in order to find the expression of the series associated to another 
region. 
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A Formulary 

In this appendix, we give some of the conventions and a few simple formulas that have been 
used to organize the calculations in this paper. 
Let us consider the function 



h{s,t) = X ^y 



nLiT'Kcfcs+^i+A) 



(130) 



where aj,bj,ej,Ck,dk and fk are real numbers and Qj and r^ are positive integers, and 
let us define 

/i(i'0)(s,i) = |-/i(s,t) and /i(°'i)(s,t) = |-/i(s,t). (131) 

/i(i'0)(s,t) = /i(s,t)A(s,t) and h^^'^\s,t) = h{s,t)B{s,t), (132) 



Then 



where 



and 



A(s, i) = - Inx + ^ QjUjipiajS + bjt + Cj) - ^ rkCkip{ckS + dkt + fk) (133) 

j=i k=l 

m p 

B{s,t) = -Iny + ^qjbjip{ajS + bjt + Bj) -'^rkdk')p{ckS + dkt + fk). (134) 

i=i fc=i 

Obviously 

m p 

A('''i){s,t) = Y,qja]+^^lj^^+'\a,s + b,t + e,) - J^rfeC^idiV^'+'^lcfcS + 4t + A) (135) 
j=i fc=i 

and 

m p 

B('='')(s, t) = ^ g,-46;+V^'+'^(«,^ + bjt + e,)-Y, rk4d[+'^^''-^'\ckS + dkt + fk). (136) 
j=i fc=i 

Now (we avoid to write the s and t dependence in the following), 

/i(i'i) = /i Tab + a(0'1)1 = /i Tab + b^^'O)] , 



since of course A^"'-*^^ = B'--*^'*^-', 



/i(0'2) = h\B^ + B(0'1) 



(137) 

(138) 
(139) 



/i(2.i) = /i [b f A2 + A(i'0)'i + 2AA(0'1) + A(1'1)1 = /i [b f A2 + A^i'O)! + 2AB(1'0) + B(2.o) 



(140) 



/i(i'2)=/i a(b2 + b(0'1) +2BB(1'°) + B(^'1) =/i A(b2 + b(0'1) +2BA(0'^)+A(0'2) 



fl41) 



etc. 



^In particle physics perturbative calculations, aj,bj,Ck and dfc are integers. 
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B Regions of convergence 



To ease the reading of this paper, we recall in this appendix some results on the theory of 
convergence of multiple (hypergeometric) series used in the text. See for instance |18| and 



19] for more details. We also give the proof in section B.4 that the polygamma functions 



ip^ ' involved in the multiple series obtained from multiple Mellin-Barnes (in this paper and, 
we think, more generally in perturbative quantum field theory) do not affect the regions of 
convergence that may be obtained, when ignoring these polygamma functions, from Horn's 



theorem (the latter being described in section B.2.1). 



B.l Single series 

In this simple case, one makes use of d'Alembert and Raabe-Duhamel's ratio tests. 

B.2 Double series 

B.2.1 Horn series 

A general statement as the d'Alembert's ratio test does not exist for arbitrary double series. 
Nevertheless, for a certain type of double series (those of Horn's type), it is possible to know 
the region of (absolute) convergence from a kind of extension of d'Alembert's ratio test due 
to Horn [TSlll^. 



A double series > > o,m,nx"^y^ is a Horn's one if the two functions 

Q(m, n) 



f{m, n) 



m=0 n=0 

Qm+i.n ^ F{m,n) 
am,n R(m, n) 



and 



l{m,n) 



a. 



m,n+l 
0"ni.'n. 



S(?7i, n) 



^142) 



are rational functions. This means that P, Q, R and S are polynomial of degree p, q, r and s. 
Defining the functions 



F{m,n) = lim f{r]m,rin) 
one may distinguish five different cases: 



and G{m,n) = lim g(r]m,r]n) . (143) 



Conditions Region of convergence 



p > r or q > s 



\x\ = and \y\ = 



p < r and q < s (|x|, \y\) G R^ 

p < r and q = s \x\ & IR+ and \y\ < ,q,q ^. 



p = r and q < s \x\ < \p(Iq\\ and \y\ G R 



p = r and q = s 



{\x\,\y\)ecnv 
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where the domains C and T) are two subsets of R^ defined as 

C = {(|x|, |y|) I < \A < ^^ and < |y| < ^^| (144) 

and 

V=\{\x\,\y\) I y{m,n)e Rl:0<\x\ < ^ or < \y\ < /^ }. (145) 

[ ' |F(m,n)| |G(m,n)| ' 



We recall that the logical disjunction 'or' in (145) implies, among others, that if < |x| < 



is true and < lyl < Trrr-^ — n is also true, then f < Ixl < tpt-^ — ^T or < \v\ < tttt-^ 

1^1 ( ^ ( TTJ. r). H ' \ II h I m, rj 1 1^1 (-Jim 



|F{m,n)| ^^^^^^^'^^^ ^ \y\ ^ |G(m,n)| ^^ <^^'=^ "^ ^^, '^^^^^^ \^" ^ l-^l ^ |F{m,n)| ^^ " ^ IWI ^ |G(m,n)| 

is true. 

The case where p = r and q' = s is the so-called Horn's theorem |19) . 

Let us conclude this section by saying that if am,n is a ratio of gamma functions of the 
type r(a?7i + bn + c) where a, b are integers and c is a constant (depending possibly on the 
dimensional regularisation parameter e), then obviously the corresponding double series will 
be of Horn's type. This is what we observe in the perturbative calculations of particle physics. 



modulo polygamma functions (see section B.4). 



B.3 Three dimensional case 

Horn's theorem may be generalized to the case of triple series or even higher order series |19) 
(but it is not as trivial as what could be thought at first sight, since it seems that wrong results 
have been given on this subject in a certain number of papers of the specialized mathematics 
litterature [20]). 

oo 

Let us consider the series > a^^^^^x^^y'^z where 

m,n,l=0 

i{m,n,tj = '-^ , g{m,n,i) = — '■ , h(m,, n,^) = — '— (146) 

0-m,n,£ (^m,n,£ (^m,n,£ 

are rational functions. We focus on the particular situation where the degree of each numerator 
is equal to the degree of its corresponding denominator. 
From the definitions 

F{m,n,i) = lim l{r]m,r]n,r]i) , G{m,n,i) = lim g{r]m,r]n,r]i) 



and 



H(m,n,£)= lim h(r/r?i, r/n,r/£) , (147) 



the region of convergence of the triple series is given by 

{\x\,\y\,\z\)eCnVinV2nV3nV4 , (148) 
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where 



e 



V, 



m,\y\,\z\ 



m,\y\,\z\ 



^\x\ < -^-, TT and ^ lyl < 7-— 

' ' |F(1,0,0)| ^'^' |G(0,1, 



and ^ 1 2:1 < 



1 



|H(0,0,1) 



\/{m,n,i) £ ]R,+ ,|x| < 



1 



\F{m,n,i)\ 



or \y\ < 



1 

\G{m,n, 



(149) 



or \z\ < 



V2 

and 
P4 



m,\y\,\z\ 



m,\y\,\z\ 



m,\y\,\z\ 



|H(m,n,^)| I' 



V(n,£) GRl,\y\ < 



y{m,£) e]R+,|x| < 



1 

|G(0,n,^)| 

1 

|F(m,0,£)| 



or \z\ < 



or \y\ < 



|H(0,n,£)|j' 
G(m, 0,^)1 f' 



y{m,n) G Ri, |x| < 



\F{m,n, 



or \z\ < 



1 



|H(m,n,0)| 



(150) 

(151) 

(152) 
(153) 
(154) 



B.4 "Almost" Horn series 

It is clear that although the integrands of the Mellin-Barnes integrals considered in this paper 
are products of gamma functions, the coefficients which appear in their (multiple) series 
representations are not purely composed of gamma functions in general. They contain also 
polygamma functions tp^ ' whose arguments contain both summation indices m and n (as 
well as e, the dimensional regularisation parameter, in the case of the box integral), implying 
that the double series are not strictly of Horn's type. 

In fact, in full generality (see Appendix Q, we observe that the coefficients of the series 
obtained from a twofold Mellin-Barnes integral whose integrand is a ratio of gamma functions, 
cannot be anything else than a ratio of gamma functions multiplied by polygamma functions 
of different orders. 

Let us call am,n this ratio of gamma functions (with arguments of the type am + bn + c 
where a and b are integers and c is not a negative integer), then 



F{x,y) = Y^ am,nx'^y'^ 



(155) 



m,n=Q 



is a Horn series, converging (absolutely) in a region TZ. 

Let ^ be a linear combination of m and n such that /i £ M, and let us define 



F{x,y) 



m,n=0 



a™,„V'^P)(a + /x)x'"y", 



(156) 



where p is a positiver integer and a is not a negative integer. We want to show that F{x, y) 
converges absolutely in the same region TZ as F{x,y). For this we follow the same reasoning 
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than the one used in Chapter 4 of |19) to prove that the region of convergence of a hypergeo- 
metric series is independent of the parameters. 
It is straightforward to prove that 

i;(p\a + fi) = 0{fi). (157) 

Then, one can find positive numberq^Co, Ci and C2 such that, for e > and /i G M, 

Co < |^(^H« + fi)\<Cifi< C2(l + e)™+" . (158) 

Thus one has 

CoFabsXx,y) < Fabs.{x,y) < C2Fabs. (x(l + e), y(l + e)) , (159) 

where the index abs. is meant for the series with absolute values. 
Now, for (x, y) E 7^, one can choose £ such that 

(x(l + e),y(l + e))G7^. (160) 

This ends the proof that F is absolutely convergent on TZ. 

This result is straightforwardly extended to multiple series of the same type with an 
arbitrary number of variables. It may also be used iteratively to prove that a multiple series 
where the general term has an arbitrary product of polygamma functions of any order do have 
the same region of convergence TZ than the series without these polygamma functions. 

Therefore one may conclude that in our cases of study, the region of convergence of a 
given multiple series may be obtained by ignoring polygamma functions in its general term, 
and by considering only the term am^n which is made of gamma functions. 

B.5 Generalized Kampe de Feriet series 

A useful class of multiple hypergeometric series, called generalized Kampe de Feriet series, is 
defined as [T9] 

00 00 ki b 

F?;:^;::::^„ (xi, . . . ,x„) = 5; • • • j; a,,_C-^. • • ^ (lei) 

A:i=0 fc„=0 "■ 



where 



n Qi 



W{aj)k,+-k,,WW{hij)k, 



"^1.-.^" = -I ^n^, • (162) 

JJ(aj)fc,+...fc„ JJJJ(A,j)fc. 

For convergence of these series one has to check that 

1 + I + nii— p — qi>Q, i = 1, . . . , n. (163) 

Moreover, the region of convergence is given either by 

n 

^|xj|^ < 1 if p> I (164) 

or by 

max{|xi|, . . . , |xn|} < 1 if p < i- (165) 



""For the case p = 0, one has to check that a + yu is such that ip{a + fj.) ^ 0. 
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